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In  experiments  conducted  on  the  NRL  Gamble  II  device, 
focused-ion  currents  have  previously  been  inferred  by  combining 
time-integrated  neutron  yield  data  with  time-resolved  diode 
electrical  characteristics  in  a model  for  ion  production.  The 
time  variation  of  ion  current  can  be  determined  in  a more 
direct  fashion  from  neutron  time-of -flight  measurements  at  — - 
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Temporal  Deuteron  Current  Determinations 
Using  Neutron  Time-of -Flight 

I.  INTRODUCTION 

The  development  of  intense  ion  beams  from  high-power 
C101ZW)  short-duration  ( 10_8s)  pulsed-plasma  sources  has  re- 
quired new  approaches  to  diagnosing  such  beams.  Measurements1 
of  intense  pulsed  proton  beams  from  the  NRL  Gamble  II  pulser 
have  been  carried  out  using  nuclear  activation2  and  ion  time- 
of-f light  with  a scintillator-photodiode  detector.  These 
techniques  have  not  been  completely  satisfactory  due  to  the 
intense  burst  of  energy  associated  with  the  proton  beam.  The 
nuclear  activation  technique  was  limited  by  the  ablation  of 
the  surface  of  the  target  being  activated.  Only  the  use  of 
attenuating  screens  eliminated  this  problem  for  proton  inten- 
sities exceeding  1018  protons  per  pulse.  The  time  variation 
of  proton  current  has  been  measured  with  a scintillator-photo- 
diode detector  in  a time-of-f light  arrangement,  but  it  is 
difficult  to  provide  a sufficiently  thin  light  shield  on  the  , 

scintillator  that  can  survive  these  ion  intensities  intact. 

Furthermore,  magnetic  or  electrostatic  deflection  is  compli- 
cated because  the  beam  is  charge  neutralized.  Both  of  these 
diagnostics  are  even  more  difficult  to  use  with  the  develop- 
ment of  more  intense  ion  beams  and  focused  beams.3  For 
focused  beams  of  deuterons,  neutron  measurements  have  been 
used  to  determine  beam  intensities  exceeding  1017  deuterons/ 
pulse.3  In  this  report,  we  present  results  of  a neutron 
diagnostic  using  the  time-of-f light  ( TOF)  technique  to  provide 
temporal  information  about  the  ion  current  for  intense  pulsed 
deuteron  beams. 

Manuscript  submitted  July  3, 1978. 
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Deuteron  beams  are  directed  onto  a deuterated  polyeth- 
ylene target  to  produce  neutrons  by  the  d(d,n)3He  reaction. 

The  neutrons  are  measured  in  time  for  different  flight  paths 
to  determine  the  time-of -emission  of  the  deuterons  relative 
to  the  diode  voltage  pulse  and  to  evaluate  the  time  variation 
of  the  deuteron  current  in  the  beam.  The  energies  of  deuter- 
ons in  the  beams  are  less  than  1 MeV  so  neutrons  produced  by 
the  12C(d,n)13N  reaction  are  separated  in  time  from  the  d-d 
neutrons.  Total  yields  of  up  to  1011  neutrons  are  measured 
on  a single  shot  with  a rhodium  activation  detector.4  For  TOF 
measurements,  fast-recovery,,  high-linearity  scintillator- 
photomultiplier  detectors  are  used.  The  scintillators  are 
6.7-cm  dia  by  5.6-cm  thick  NE-111*  which  are  mounted  on  XP- 
2020  photomultiplier  tubest.  The  time  resolution  ( FWHM)  of 
these  detectors  to  individual  neutron  interactions  is  about 
5 ns. 

II.  DEUTERON  EMISSION  TIME 

To  determine  the  neutron  emission  time,  two  TOF  detec- 
tors at  15°  to  the  axis  of  the  diode  with  different  flight 
paths  of  1.91  m and  7.42  m are  used.  Each  detector  is  enclosed 
in  a lead  shield  with  2.5  cm  of  lead  in  the  flight  path  to 
reduce  the  saturation  level  of  the  x-ray  flash  in  the  scintil- 
lator. In  addition,  a neutral  density  filter  is  used  between 
the  scintillator  and  photomultiplier  to  prevent  saturation  of 
the  neutron  response  in  the  photomultiplier.  The  responses  of 
these  detectors  on  a particular  shot  are  shown  in  Fig.  1. 

The  detectors  recover  from  the  initial  x-ray  flash  before  the 
neutron  responses  are  measured.  The  neutron  response  corre- 
sponding to  the  longer  flight  path  is  delayed  in  time  as 
expected. 

* 

Available  from  Nuclear  Enterprises,  Ltd.,  San  Carlos,  CA,  94070 
^ Available  from  Amperex  Electronic  Corp.,  Hicksville,  NY,  11802 
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The  energy  and  emission  time  of  the  most  energetic 
neutrons  are  deduced  from  the  time-of-f light  traces  in  Fig.  1. 
The  most  energetic  neutrons  have  a flight  time  corresponding 
to  the  leading  edge  of  the  neutron  response.  We  use  time  on 
the  leading  edge  corresponding  to  one-half  of  the  signal  max- 
imum. It  is  assumed  that  no  high  energy  neutrons  are  emitted 
late  in  the  pulse  and  reach  the  detector  before  lower  energy 
neutrons  emitted  earlier  in  time.  This  assumption  will  be 
checked  later.  The  neutron  time  for  each  detector  is  measured 
relative  to  the  peak  of  the  x-ray  signal  and  corrected  for  x- 
ray  flight  time.  These  flight  times  correspond  to  a maximum 
neutron  energy  of  3.41  -t  0.15  MeV.  The  corresponding  deuteron 
energy  is  0.45  -t  0.10  MeV  based  on  the  d(d,n)3He  reaction  at 
15° . The  uncertainties  in  these  values  are  a result  of  uncer- 
tainties in  time  and  distance  measurements.  The  neutron 
flight  time  versus  distance  is  extrapolated  to  zero  distance 
to  determine  the  emission  time  for  the  most  energetic  neutrons. 

The  emission  time  and  deuteron  energy  are  compared  with 
the  time  dependent  accelerating  voltage  in  Fig.  1.  The  ion 
energy  trace  in  this  figure  corresponds  to  the  measured  diode 
voltage  corrected  for  the  drift  tube  inductance  and  reduced  by 
the  energy  loss  of  the  deuterons  in  traversing  a thin  foil 
used  for  current  neutralization.3  The  maximum  deuteron  energy 
for  this  trace  is  0.46  MeV  in  good  agreement  with  the  energy 
determination  by  TOF.  The  neutron  emission  time  is  determined 
relative  to  this  ion  energy  trace  through  a fiducial  timing 
pulse.  This  corresponds  to  a time  of  87  ns  on  the  ion  energy 
trace  in  Fig.  1.  This  time  must  be  reduced  by  the  flight  time 
of  deuterons  from  the  diode  to  the  target  ( 20  ns)  to  give  the 
emission  time  of  deuterons  from  the  diode  of  67  ± 7 ns  as 
shown  in  Fig.  1.  The  uncertainty  in  this  value  results  from 
propagating  uncertainties  in  time  and  distance  measurements. 

The  deuteron  emission  time  corresponds  to  the  maximum  deuteron 
energy  within  experimental  error.  Also  there  is  no  evidence 
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in  the  ion  energy  trace  of  more  energetic  deuterons  emitted 
late  in  time  which  is  consistent  with  the  previous  assumption 
that  the  most  energetic  neutrons  are  not  emitted  late  in  the 
pulse . 

III.  DEUTERON  CURRENT 

The  time  spread  of  the  neutron  response  in  a TOF  detec- 
tor is  composed  of  three  contributions:  1)  the  detector  re- 
sponse, 2^  the  energy  spread  of  the  neutrons,  and,  3)  the 
duration  of  emission  of  the  neutrons.  The  time  spread  of  the 
responses  in  Fig.  1 is  about  50  ns  which  is  much  larger  than 
the  detector  response  time.  Since  the  width  of  the  neutron 
response  does  not  increase  much  with  flight  path,  the  time 
spread  is  determined  in  large  part  by  the  duration  of  neutron 
emission  and  not  by  the  energy  spread  of  the  neutrons.  The 
duration  of  neutron  emission  then  provides  a direct  measure 
of  the  duration  of  the  deuteron  current. 

1 . Relation  of  Deuteron  Current  and  TOF  Yield 

The  time  variation  of  the  neutron  yit  ( tn)  is 

related  to  the  time  variation  of  the  deuteron  current  Nd( t) 
in  the  following  way: 

Yn(tn)  - Nd(t)  f d a/e dE  ( 1) 

J o 

• 

where  tn  - t + «,n/vn  + <-d/vd 

ln  - neutron  flight  path 

vn  - neutron  velocity 

td  - deuteron  flight  path 

Ed  “ l/2mdvd2  • incident  deuteron  energy 

a - cross  section  for  the  d(d,n)3He 
reaction 
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e = stopping  cross  section  for  deuterons 
in  polyethylene 

The  neutron  yield  is  related  to  the  deuteron  current  through 
the  reaction  cross  section  integrated  over  the  target  thick- 
ness. It  is  assumed  that  the  target  thickness  is  larger  than 
the  range  of  the  deuterons.  The  time  of  arrival  of  a neutron 
at  the  TOF  detector  ( t ^ is  different  from  the  time  that  the 

deuteron  was  emitted  from  the  diode  by  the  flight  time  of  the 
deuteron  to  the  target  ( -C-^/v^  and  the  neutron  flight  time 

(l  /v  ) . The  slowing-down  time  of  the  deuterons  in  the  target 
n n 

is  negligible  compared  to  the  neutron  flight  time.  It  would 
appear  that  a measurement  of  the  time  variation  of  the  neutron 
yield  can  be  used  to  determine  the  deuteron  current  in  time  as 
long  as  the  overtaking  of  slower  neutrons  by  faster  neutrons, 
i.e.  orbit  crossing,  is  not  too  severe. 

In  practice,  the  problem  is  complicated  by  two  effects. 
First,  the  deuterons  lose  energy  in  the  target  so  the  outgoing 
neutrons  have  a velocity  spread  and  hence  a time  spread. 

Since  the  cross  section  changes  with  deuteron  energy,  the 
neutron  yield  will  change  with  the  neutron  energy  or  flight 
time.  Second,  for  a configuration  in  which  the  ion  beam  is 
focused,  there  is  a spread  in  angle  for  the  deuterons  incident 
on  the  target  which  produces  a spread  in  angle  for  the  out- 
going neutrons.  In  our  experiments,  the  half-angle  for  the 
focusing  cone  of  deuterons  on  the  target  is  about  35° . Since 
the  cross  section  and  neutron  energy  depend  on  the  neutron 
emission  angle,  both  a time  spread  of  the  outgoing  neutrons 
and  a variation  in  neutron  yield  arise  from  this  angular 
spread. 

The  dependence  of  neutron  energy  on  deuteron  energy  for 
different  neutron  emission  angles  is  shown  in  Fig.  2.  At  small 
angles  the  neutron  energy  changes  rapidly  with  deuteron  energy, 
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but  near  90°  the  neutron  energy  is  less  sensitive  to  the 
deuteron  energy.  If  the  deuteron  beam  is  unidirectional  on 
the  target,  as  in  an  unfocused  diode  geometry,  then  the 
neutron  time  spread  can  be  minimized  by  locating  the  detector 
near  90° . 

The  dependence  of  the  cross  section  on  deuteron  energy 
and  neutron  emission  angle  is  shown  in  Fig.  3.  At  all  angles 
the  cross  section  increases  with  deuteron  energy.  Conse- 
quently, the  neutron  response  is  largest  for  the  most  ener- 
getic neutrons  and  decreases  in  intensity  toward  longer  time 
due  to  the  deuteron  energy  loss  in  the  target.  Near  90°,  the 
cross  section  is  relatively  insensitive  to  neutron  angle  so 
that  in  a focused  geometry  with  angular  spread,  all  angles 
will  contribute  nearly  equally  to  the  neutron  response. 

In  our  experiment,  not  only  do  we  have  the  energy  and 
angle  variations  discussed  previously,  but  also  the  deuteron 
energy  is  time  dependent  due  to  the  time  variation  of  the 
diode  voltage.  For  a neutron  detector  located  at  an  angle  0^ 
relative  to  the  deuteron  beam  direction,  Eq.  1 can  be  re- 
written to  include  these  effects. 

r E _( t) 

Yn(tn’V  “ Nd(t)  J0  a(E,0n)/e(E)dE  (2) 

l t . 

where  t - t + — + — . 

n 

n d 

The  neutron  velocity  is  determined  by  the  kinematics  of  the 
d(d,n)3He  reaction  for  given  values  of  Ed  and  0n>  A straight- 
forward inversion  of  Eq.  2 for  deuteron  current  is  not  possible. 
Instead,  the  right  side  of  Eq.  2 is  calculated  and  compared 
with  the  measured  neutron  response.  For  this  purpose,  the 
corrected  diode  voltage  trace  is  used  for  Ed( t)  and  the  deu- 
teron current  predicted  by  a simple  analytic  model1  is  used 
for  Nd( t)  . 
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Values  of  the  reaction  cross  section  (a)  and  the  stop- 
ping cross  section  (e)  are  extracted  from  the  literature. 
Center  of  mass  differential  cross  sections  for  the  d(d,n)3He 
reaction  from  97  keV  to  350  keV  (Ref.  5)  and  from  300  keV  to 
700  keV  (Ref.  6)  were  transformed  to  the  laboratory  system  and 
fitted  with  smooth  curves  such  as  given  in  Fig.  3.  These 
cross  sections  are  based  on  measuring  the  outgoing  3He  par- 
ticles. At  higher  energies,  cross  sections  given  by  Brolley 
and  Fowler7  are  used.  These  cross  sections,  which  are  based 
on  neutron  measurements,  are  consistent  with  values  deduced 
from  Ref.  6 below  700  keV  except  at  the  smallest  angles 
(0°  and  10°).  For  these  angles  the  values  of  Brolley  and 
Fowler  are  about  13#  higher  than  values  deduced  from  Ref.  6 
and  were  normalized  to  the  values  from  Ref.  6.  Stopping 
cross  sections  for  deuterons  in  polyethylene  were  calculated 
from  stopping  cross  sections  for  protons  in  hydrogen  and 
carbon  taken  from  the  recent  compilation  of  Anderson  and 
Ziegler.8  Values  of  a and  e which  were  used  in  evaluating 
neutron  yields  with  Eq.  2 are  listed  i.^able  1. 

2 . Measurements  and  Results 

Neutron  TOF  measurements  have  been  made  with  detectors 
located  at  15°  and  85°  to  the  incident  deuteron  beam  direction 
for  flight  paths  of  1.91  m and  2.36  m respectively.  Short 
flight  paths  are  used  to  minimize  the  neutron  time  spread  from 
energy  variations.  However,  these  flight  paths  are  suffi- 
ciently long  for  the  detectors  to  recover  from  the  x-ray  flash 
before  recording  neutrons.  Results  for  a planar  diode 
geometry  without  focusing  are  shown  in  Fig.  4.  Here  measured 
TOF  responses  are  compared  with  yields  calculated  according 
to  Eq.  2.  The  deuteron  energy  Ed( t)  and  the  deuteron  current 

Nd( t)  for  these  calculations  are  given  in  Fig.  4.  Neutron 

yields  are  calculated  in  2 nsec  time  increments  of  Ed( t)  and 

Nd( t)  and  then  summed  to  give  the  curves  which  are  plotted 


as  open  points.  The  calculations  have  been  normalized  to  the 
peaks  of  the  measured  responses. 


To  assess  the  importance  of  the  duration  of  the  current 
and  energy  pulses  on  the  neutron  response,  a calculation  was 
carried  out  using  a monoenergetic  deuteron  pulse  of  constant 
current  with  a duration  of  only  2 nsec.  Neutron  responses 
calculated  in  this  case  are  shown  at  the  bottom  of  the  figure. 
Here  the  temporal  spread  is  larger  at  15°  than  at  85°  as 
expected  from  the  cross  section  and  kinematic  variations  dis- 
cussed previously.  However,  at  both  angles  these  contribu- 
tions are  much  less  than  the  contribution  from  the  duration 
of  ion  emission  as  given  by  the  ion  current  and  energy  traces. 

Results  for  the  focused  diode  geometry  are  given  in  Fig. 
5.  In  this  geometry,  the  neutron  angles  range  from  0°  to  50° 
for  the  15°  detector  and  from  50°  to  120°  for  the  85°  detector. 
Neutron  yields  are  calculated  in  10°  increments  over  these 
ranges  according  to  Eq.  2 and  summed  to  give  the  curves  which 
are  plotted  as  open  points  in  Fig.  5.  The  deuteron  energy 
Ed( t)  and  deuteron  current  Nd( t)  used  for  these  calculations 
are  given  in  Fig.  5.  Again  the  calculations  have  been  normal- 
ized to  the  measured  responses.  Responses  calculated  for  a 
monoenergetic  deuteron  pulse  of  2 nsec  duration  are  shown  at 
the  bottom  of  the  figure.  The  larger  temporal  spreads  of 
these  responses  than  those  in  Fig.  4 indicate  that  the  angular 
spread  is  limiting  the  temporal  resolution  significantly.  Now 
the  temporal  spread  is  larger  at  85°  than  at  15°  as  expected 
from  the  angular  variations  of  the  cross  sections  and  kine- 
matics discussed  previously.  Even  with  these  temporal  reso- 
lutions, the  duration  of  the  neutron  T0F  traces  results  pri- 
marily from  the  duration  of  ion  emission  as  given  by  the  ion- 
current  and  energy  traces. 

In  making  the  comparison  between  calculated  yields  and 
measured  T0F  responses,  it  is  assumed  that  the  neutron 
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response  is  directly  proportional  to  the  neutron  yield  given 
by  Eq.  2.  The  neutron  response  Rn  is  related  to  the  neutron 

yield  Yr  by  the  following  expression. 

Rn(t)  = n e P Yn  ( t) 

where  Q = solid  angle  of  the  neutron 

detector 

ft  = efficiency  of  the  neutron 
detector 

P = pulse-height  response  of  the 
detector  to  neutrons. 

The  quantities  ft  and  P are  complicated  energy-dependent 
properties  of  the  neutron  detector.  The  primary  interaction 
process  in  the  scintillator  is  n-p  elastic  scattering.  In 
this  case,  the  efficiency  is  proportional  to  the  total  n-p 
cross  section  and  the  pulse-height  response  depends  on  the 
proton-recoil  energy,  which  is  in  fact  a pulse-height  distri- 
bution due  to  the  variation  in  scattering  angle  of  the  recoil 
protons.  An  average  over  this  distribution  is  made  in  the 
experiment  because  ~ 10°  neutrons  interact  in  each  detector 
in  a single  shot.  Then  the  quantity  P corresponds  to  a pulse- 
height  response  averaged  over  the  recoil-proton  pulse-height 
distribution.  Since  this  pulse-height  response  increases  with 
neutron  energy  and  the  n-p  cross  section  decreases  with 
neutron  energy  in  the  energy  range  of  interest,  these  varia- 
tions tend  to  cancel  each  other.  In  the  flat  diode  geometry, 
it  is  estimated  that  the  net  variation  in  the  product  of  ft 
and  P is  £ 7#  for  Ed  ranging  from  500  to  700  keV.  In  the 

focused  diode  geometry  with  a larger  angular  spread  of  the 
neutrons,  the  net  variation  in  the  product  of  ft  and  P is  15# 
for  the  15°  detector  and  30#  for  the  85°  detector  for  the  same 
range  of  deuteron  energies.  Neutron  interactions  with  carbon 


in  the  scintillator  have  been  ignored  in  making  these 
estimates . 


IV.  CONCLUSIONS 

Within  experimental  uncertainty,  the  time  of  emission 
of  the  most  energetic  neutrons  corresponds  to  the  production 
of  the  most  energetic  deuterons  in  the  diode.  In  both  the 
focused  and  unfocused  geometries,  the  duration  of  the  neutron 
response  measured  at  either  15°  or  85°  to  the  anode  to  cath- 
ode direction  results  primarily  from  the  duration  of  deuteron 
emission  in  the  diode.  In  the  focused  geometry,  the  optimum 
time  resolution  is  achieved  at  a small  angle,  but  in  the 
unfocused  geometry,  the  time  resolution  is  optimum  at  an 
angle  near  90°.  Finally,  the  shapes  of  the  neutron  yields 
calculated  from  an  experimentally  determined  ion  energy 
trace  and  a model  dependent  ion  current  trace  are  in  reason- 
able agreement  with  the  measured  responses  for  both  focused 
and  unfocused  geometries  at  both  angles.  This  last  result 
is  especially  important  because  time  dependent  ion  currents 
are  difficult  to  measure  directly  for  the  intense  beams  of 
interest.  The  agreement  between  the  measured  and  calculated 
neutron  yields  supports  the  validity  of  the  ion-production 
model  used  to  unfold  the  ion  current  from  measured  electri- 
cal characteristics  of  the  diode. 
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MAXIMUM  DEUTERON  ENERGY  • (0.45  ±0.1 ) MeV 


> 


TIME  OF  EMISSION  OF  MAXIMUM  ENERGY  DEUTERONS  - (67  ±7)n««c 

Fig.  1 — Neutron  TOF  measurements  for  two  detector*  at  15°  to  the  diode  axis 
and  at  flight  path*  at  1.91  m and  7.42  m respectively.  The  maximum  deuteron 
energy  corresponds  to  the  moat  energetic  neutron*  which  are  detected  at  the  leading 
edge  of  the  TOF  trace*.  On  the  bottom  of  thia  figure,  the  emission  time  of  these 
deuterona  from  the  diode  is  compared  with  the  ion-energy  trace. 
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NEUTRON  ENERGY  (MeV) 


Fig.  2 — Dependence  of  outgoing  neutron  energy  on  incident  deuteron 
energy  for  the  d(d,n)3He  reaction  for  neutron  emission  angles  ranging 
from  0°  to  120°. 
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Fig.  3 — Dependence  of  the  d(d,n)3He  reaction  cross  section  on 
incident  deuteron  energy  for  neutron-emission  angles  ranging  from 
0°  to  120°. 
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Fig.  4 — Remit*  of  neutron  TOF  measurements  at  15°  and  85°  for  a flat  (unfocused) 
diode  geometry.  The  diode  geometry  is  shown  at  the  top  left:  A * Anode,  K * 
Cathode,  T - Thin  foil.  The  ion-energy  and  current  traces  for  the  shot  are  shown  at 
the  top  right.  The  measured  neutron  TOF  traces  (solid  lines)  are  compared  with  cal- 
culated yields  (open  points).  Calculations  for  a monoenergetic  deuteron  pulse  of 
2 nsec  duration  are  given  at  the  bottom  of  the  figure. 
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FOCUSED  GEOMETRY 


Fig.  B — Results  of  neutron  TOF  measurements  at  15°  and  85°  for  a focused-diode 
geometry.  The  diode  geometry  is  shown  at  the  top  left:  A * Anode,  K “ Cathode, 
T * Thin  foil.  The  ion-energy  and  current  traces  for  this  shot  are  shown  at  the  top 
right.  Measured  neutron  TOF  traces  (solid  lines)  are  compared  with  calculated  yields 
(open  points).  Calculations  for  a monoenergetic  deuteron  pulse  of  2 nsec  duration 
are  given  at  the  bottom  of  the  figure. 
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